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The geometric phases of cychc evolutions for mixed states are discussed in the 
framework of unitary evolution. A canonical one- form is defined whose line integral 
, gives the geometric phase which is gauge invariant. It reduces to the Aharonov and 

' Anandan phase in the pure state case. Our definition is consistent with the phase 

m : 

■ shift in the proposed experiment [Phys. Rev. Lett. 85, 2845 (2000)] for a cyclic 
\ evolution if the unitary transformation satisfies the parallel transport condition. A 

comprehensive geometric interpretation is also given. It shows that the geometric 

■ phases for mixed states share the same geometric sense with the pure states. 
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When a pure quantal state undergoes a cyclic evolution the system returns to its original 
state but may acquire a nontrivial phase factor of purely geometric origin. This was first 
discovered^by Berry [Ij in the adiabatic context, and generalized to non-Abelian by Wilczek 
and Zee A nice interpretation was given by Simon jsl in terms of a natural Hermitian 
connection, as the parallel transport holonomy in a Hermitian h.ne bundle. Extension to the 
nonadiabatic cyclic case was given by Aharonov and Anandan [4]. Based on Pancharatnam's 
earlier work on interference of light Isf , this concept was generalized to noncyclic evolutions 

no 

and nonunjtary evolutions p, 13| . The geometric phases for entangled states have also been 



discussed ISl. Applications of the geometric phase have been found in molecular dynamics 

Q n 

[9|], response function of the many-body system [10|], and geometric quantum computation 



In all these developments the geometric phases have been discussed only for pure states. 
However, in some applications we are interested in mixed state cases 

HQ. 

Uhlmann was probably the first to address the issue of mixed state in the context of 
purification, but as a purely mathematical problem Recently, Sjoqvist et al. 



gave a 



new formalism of the geometric phase for mixed state in the experiment context of quantum 
interferometry under parallel transport condition. It has been pointed out 1J| that the 
latter geometric phase can be undefined at nodal points in the parameter space where the 
interference visibility vanishes. Anyway, the geometric phases for mixed states proposed in 
Refs. and are generically different in the unitary case and match only under very 
special conditions such as in terms of pure states ,15]. 

In this paper, we give a definition of geometric phase for a cyclic evolution of mixed 
quantal state in the dynamical context of quantum system. The reasons of employing the 
cyclic evolution are two: (i) the cyclic evolution of a physical system is of the most interest in 
physics both experimentally and theoretically (ii) the phase shift in a cyclic evolution should 
be definite. Firstly, we give a straightforward generalization of Aharonov and Anandan (A- 
A) phase for the global cyclic evolution where the total phase is explicit. Though this case 
seems a trivial extension of (A-A) phase, it contains the essence of geometric phase of mixed 
state. Then, we give the discussion of the general case based on a definition of the total 
phase. This geometric phase reduces to the (A-A) phase P], the standard geometric phase 
for pure state undergoing a cyclic evolution. We find that if the evolution satisfies the 
parallel transport condition the geometric phase is consistent with the result in Ref. jl^ . 
Moreover, we give the geometric meaning of the geometric phases of mixed states which 
share the same sense with the pure states for the first time. 

Supposing a quantum system with the Hamiltonian H{t), the density operator p{t) of 
this system will undergo the following evolution 

p(t) = f/(t)p(0)f/+(t), (1) 

where U{t) = Te^^-fo is a unitary transformation, here T is the chronological oper- 

ator. If f/(r) and p(0) are commutative: [[/(r),p(0)] = 0, i.e., p(r) = p(0), we say this state 
undergoes a cyclic evolution with period r. Furthermore, if U{t) = e^'^I, this evolution can 
be called the global cyclic evolution since for any p(0) we have p(r) = p(0). 

At first, we study the case of the global cyclic evolution. Now define U{t) = e~*'^*^*^f/(t) 
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such that U{t) = I. We define the geometric phase for such state during the global cyclic 
evolution as 



Tr 



p(0)t7+(t) 



dU{t) 
dt 



dt. 



Using the transformation between U{t) and U{t), one can have 



and 



-I I dtTr 



dU{t) 
dt 



(2) 



(3) 



(4) 



-I £ dtTr [pit)Hit)]. 



Obviously 0^ is just the dynamical phase during the cyclic evolution. 

We can prove that if p{t) is the density operator of a pure state the geometric phase defined 
by Eq. Q is just the Aharonov and Anandan phase 4]. Assuming p(0) = {ipiO)) ("0(0)1 , 
\^{t)) = U{t) 10(0)). Let \^{t)) = e-*'^^*) |0(t)) , then we have \^{t)) = U{t) |(^(0)) . So, the 
Eq. (121) can be written as (pg = i Jq {(f{t) \ 0(t)) dt which is just the result of Ref. j^. 

An initial state can always be diagonalized, namely, 

p{0)=J2w,\k){k\, (5) 

k 

where \k) are bases for the system and Wk are classical probability to find a member of the 
ensemble in the corresponding state. For the global cyclic evolution we have U{t) \k) = e 



i(f> 



\k) , then the A-A phase of \k) is 



where 0^ = -i £ dtTr \\k) {k\ U+{t)U{t) 



Substituting into (jlj), and from ^ we obtain 



(6) 



k k 

So, for the global cyclic evolution the geometric phases of mixed states have explicit mean- 
ings: the geometric phases of a mixed state is the weighted average of the geometric phases 
of the constitute pure states. 

Example I. -Suppose that a qubit (a spin-| particle) with a magnetic moment is in a 
homogenous magnetic field B along the z axis. Then the Hamiltonian in the rest frame is 
H = —pBaz- Suppose the initial state is 

1, 



P(0) = -[I + r{sinea^ + cosOa,)], 



(7) 
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where r is a constant and < r < 1. So we have p{t) = U{t)p{0)U~^{t) with 

U{t) = exp{ifiBtaJh). (8) 

This unitary evolution is periodic with period r = nfi/pB, i.e., p(r) = p(0). It is easy to see 
that U{t) = exp(z7r)/. Let U{t) = e-^^-^*/^f/(t), then U{t) = U{0). From Eq. (0), and after 
some elaboration we can obtain the geometric phase 



^9 



TX 



{l-rcosO). (9) 



Obviously, if r = 1, we get (j)g = 7r(l — cos 6') which is just the A-A phase 

On the other hand, we can have two pure states pi(0) = |[/ + (sin^cr^. + cos^cr^)] and 
P2(0) = |[/ + (sin(7r + 9)ax + cos(7r + 9)az)], which can construct a set of orthonormal bases. 
From 0, we have p(0) = ^pi(O) + ^P2(0). Obviously, this is a diagonal representation 
of the initial state p(0). Then from Eq.®, we can obtain (pg = ii^7r(l — cos^^) + ^^7r(l — 
cos(7r + 9)) = 7r(l — r cos 6*). 

The above discussion of the general cyclic evolution seems a trivial extension of (A-A) 
phase, but it contains the essence of geometric phase of mixed state. 

For the general case of a cyclic evolution, the density matrix and transformation satisfy 
[f/(r), p(0)] = 0. We can not find the total phase explicitly from this condition. To factor out 
ihe total phase, we use the Pancharatnam's brilliant idea, i.e., the Pancharatnam connection 



s 



We define the total phase of the mixed state during a cyclic evolution with the initial 
state p(0) and the unitary transformation U {t) as 

= argTr[p(0)f/(r)]. (10) 

Let U{t) = e~^'^^^^U{t) such that 0(r) = 0. Based on this definition, the geometric phase 
of the cyclic evolution can be also defined by Eq. (jSJ. Obviously, the geometric phase for 
a cyclic evolution takes the same form as in Eq. (jH)). Indeed Eq. Q defines a canonical 
one-form in the parameter space: 

p = iTr [p{0)U+{t)dU{t)] . (11) 

It is not difficult to prove that /? is a real number. The geometric phase can be obtained by 
its line integral, i.e., 

V. (12) 



The equivalent of the above formula for the pure states case is well-known j7|, |lf 

The geometric phase defined above is manifestly gauge invariant: it does not depend on 
the dynamics, but it depends only on the geometry of the close unitary path given by the 
unitary transformation U{t). Assuming a is a dynamic parameter of this system, the nature 
of the cyclic evolution requires a(0) = ^(t) Jl^. Under the transformation U(t)' = e*'^*^"^?7(t), 
we can have (3' = (3 — d6. It is easy to prove that § f3' = § [3 since (5(a(r)) = 5(a(0)). Indeed 
the quantity /3 = iTr[p{Q)U^{t)dU{t)\ can be regarded as a gauge potential on the space of 
density operators pertaining to the system. 

Example II. -Consider a spin-| particle is initially in the state 



p(0) = -(/ + ra. 



(13) 



B(t) 



where r < 1 is a constant. Suppose this particle is in a magnetic field B()f:) with 

— {ujh/fi)ey, < t < ti 

-{iuh/fi)ez, h <t <t2 , 

-{uh/ fi){smipex — cosipCy), t2 <t < r 

in which u; is a constant and ^1 = ^)^2 = and r = So, the unitary transformation 



(14) 



IS 



U{t) 



^—iiOtCTy 



-iu){t — t2){sirv (^(Jx—COS ^CTy) 



_■ ip ■ 

^ g ^ "2" g 'i~0''U 



0<t<ti 

tl<t<t2 

h<t<T 



(15) 



Then, U{t) 



-i I (sin ificr^ -cos ipcTy ) -i^a-^ -i 



e-h^v. We can prove that f/(r)p(0)f/+(r) = p(0). 



Under this unitary transformation, the state undergos a cyclic evolution with a closed path 
of the corresponding Bloch vector as showing in Fig.l, where the vectors at points B and 
C are = (r sin^, 0,r cos6') and rc = {r sin 6 cos ip , r sin 6 simp , r cos 6) respectively. From 
Eq. (jlUj) and (j3]), we have the total phase of this cyclic evolution: (p = — arctan rtan^ , 
and the dynamic phase: (pd = — frcos^. Then, the geometric phase of this cyclic evolution 
is 



— arctan 



r tan ■ 



+ —r cos 6. 



(16) 



We know that if 6^ = | the closed path on the Bloch sphere is geodesic [17|, i.e., the dynamical 
phase is zero. At this time the geometric phase is — arctan rtan^ , which has also been 
pointed out in Ref. Q] and verified by the recent experimental observations 3]. 
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FIG. 1: Path with cydic evolution of the Bloch vector. The sohd hne represents a geodesic path 
of the unitary pure state case and defines a spherical triangle enclosing the solid angle 7r/2. The 
dashed line represents the closed path of the mixed state with the unitary evolution defined by 
(|T5|l . where the initial state is given in Eq. (|13j) . 

If (f = 27r, i.e., the points B and C in Fig. 1 are the same. Obviously, the geometric 
phase for such a path should be equal to the case in Example I except for the direction of 
magnetic field reversed. From (fT^ . we have (pg = 7r(l + rcos^) which is consistent with Eq. 
®. 

From our definition of the geometric phase, we can get a theorem for a composite quantum 
system. 

Theorem: Let p^^ be a density operator of a composite quantum system of A and 
B. If the system evolves under the unitary transformation U{t) = with = 

Te~*-^o ^"^^ and [p^^, U{r)] = 0, then the geometric phase of such a composite system 
equals to one of the subsystem B, i.e., (j)'^^ = (pg. 

The density operator of such system can be expressed by 

+/^®r^- A^ + A,Af ® Af], (17) 
where and A'"^ are the orders of density matrices for each subsystems, A"^ = (A^; i = 
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1, 2, ■ ■ ■ , iV| - 1) and A-^ = (Af ; i = 1, 2, ■ ■ ■ , iV| - 1) are the generators of SU{Na) and 
SU{Nb) respectively, = {rf; z = 1, 2, ■ ■ ■ , iV^ - 1) and = (rf ; z = 1, 2, ■ ■ ■ , - 1) are 
two vectors, and Pij are (A^^ — 1)(A''^ — 1) real numbers. Because Tr{Xf'^) = 0, one can 
have (p^^ = aigTr[p^^U{T)] = argTr[p^f/^(r)] where is reduced density operator for 
B, i.e., = r^-A^]. From 0, we can also obtain (pj^ = -i j Tr[p^^U+{t)U{t)\ = 

-i JTr[p^U^^ {t)U^{t)]. So the geometric phase for p^^ is = - (f)^^ = (f)^ - 0f , 
where 0^ = arg Tr[p-^[/^(r)] and 0f = —i J Tr[p^U^^ (t)if^ (t)]. The geometric phase of 
the subsystem i? is 0^ = 0^ — 0f = 4>f^- 

We know that for a pure state of an entangled composite quantum system p^^, the 
reduced density operators p^'^ are mixed states. So this theorem may have some important 
applications, at least to observe the geometric phase for mixed state, since the geometric 
phase for a pure state can be acquired by methods known before. 

Recently, Sjoqvist et al. have also proposed geometric phases for mixed states under 
the parallel transport condition If the unitary transformation U{t) satisfies a parallel 
condition we can prove the geometric phase given by Eq. (j21) is consistent with what was 



12j | . The parallel transport condition for a mixed state undergoing unitary 



obtained in Ref. 
evolution is 

Tr[p{t)U+{t)U{t)] = 0. (18) 



12| | the geometric phase of the cyclic 



Under this condition, from the definition in Ref. 
evolution can be expressed as 0^ = argTr[p(0)f/(r)]. 

Indeed, for the unitary operator U{t) we can always obtain another unitary operator 
by a f/(l) transformation, U'{t) = e^«(*)f/(t), so that Tr[p{t)U'+ {t)U' (t)] = 0. From this 
condition, it is easy to prove that 

iit)=iTr[piO)U+it)Uit)]. (19) 

So ^{t) = From ^ and ©, 0^ = argTr[p(0)f/(r)] + ^(r) = argTr[p(0)[/'(r)]. 

Obviously, such a parallel transport U'{t) is unique for an initial state p(0) and the unitary 
operator U{t). So, if U{t) satisfies the parallel transport condition, the dynamical phase 
■„ E,. O vanishes .deuticaHy. so ^ .e., unde. the parallel condition ou. 

definition of the geometric phase is consistent with the definition in Ref. |12l |. 

However, for an initial state p(0) and the unitary operator U{t), one can construct al- 
ternative parallel transport by a unitary transformation, except for a U{1) transformation, 
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i.e., the parallel transport corresponding the initial state and the unitary operator U{t) 
is not unique We must emphasize that only for the parallel transport obtained 

by the U{1) transformation (which is unique), U'{t) = e*^*^*^f/(t) determined by Eq. ^T^ . 
the geometric phase can be expressed as (pg = argTr[p(0)t/'(r)]. As an example, we con- 
sider the case of Example I. For the state p(0) given in 0, one can have two pure states 
p± = |[J ± {sinOax + cos9az)] such that p(0) = + ^-f-P-- For the unitary opera- 

tor U(t) = exp{iiiBtaJfi), by unitary transformation e~*[^'^*(™''^)(''™^'^^+'=°''^'^^)/^], we can 
construct U" = e»(A'^*'^-/^)e-*[^-^*(™'^^)('^''^^'"-+™"'''^-)/^l which satisfies Tr[p±{U")+U"{t)] = 0. 
Obviously, this is a parallel transport of p(t) = U"{t)p{0)U"^ (t) = U (t) p{0)U+ (t) . It is 
easy to obtain argTr[p(0)t/"(r)] = vr — arctan[r tan(7r cos 6*)] which does not equal to the 
geometric phase given by Eq. ©. 

In fact, the set of the density operators constructs a base space Ai (the subspace corre- 
sponding to pure states is isomorphic to the projective Hilbert space). The unitary trans- 
formations U just define a line bundle JF on the base space with the projection map 11 as 
n~^(p) = {p : p ^ c = Tr{Up)/\Tr{Up)\, c G U{1)}. Then the cychc evolution discussed 
in the present paper is just define a closed curve C : [0,t] ^ A4 on the base space, but on 
the bundle space the path is not closed with the final point lifted by a U{1) factor c(r). On 
the other hand, using U{1) transformations, we can always obtain a unitary transformation 
U'{t) determined by Eq. (fTI?|) which satisfies the parallel transport condition for a given 
closed path C. Then, the geometric phase defined in Eq. can be expressed by the phase 
factor of the parallel lift : (pg = argTr[p(r)f/'(r)]. Under this consideration the geometric 
phases of mixed states share the same geometric sense with the pure states for the cyclic 
evolution. 

In summary, we give a definition of the geometric phase for a mixed state during a cyclic 



12 1 when 



evolution. Our definition is consistent with the prescription proposed in Ref. 
the evolution satisfies parallel transport condition. We first give the geometric meaning 
of the geometric phases for mixed states which share the same sense with the pure states. 
Although our definition is only for cyclic evolution, the discussion can be straightforwardly 
generalized to noncyclic evolution except for a gauge condition [2^ (for the pure states see 
in 0]) which is automatically satisfied in cyclic evolution. 
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